Collective Motion of Vibrated Polar Disks 
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We experimentally study a monolayer of vibrated disks with a built-in polar asymmetry which 
enables them to move quasi-balistically on a large persistence length. Alignment occurs during 
collisions as a result of self-propulsion and hard core repulsion. Varying the amplitude of the 
vibration, we observe the onset of large-scale collective motion and the existence of giant number 
fluctuations with a scaling exponent in agreement with the predicted theoretical value. 
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The reeent surge of theoretical/numerical activity 
about the collective properties of interacting self- 
propelled particles has produced some striking results, 
even in the simplest situations where local alignment, 
the only interaction, is competing with some noise: for 
instance, true long-range order may arise in two di- 
mensions, yielding collectively-moving ordered phases en- 
dowed with generic long-range correlations and anoma- 
lous "giant" number fluctuations [iHl] . Despite the ubiq- 
uity of the collective motion, observed at all scales in 
more or less complex situations ranging from the coop- 
erative action of molecular motors |5| , the collective dis- 
placement of cells 0, to the behavior of large animal, hu- 
man, or robot groups Q , there is, as of now, a lack of well- 
controlled experiments to which this theoretical progress 
can be seriously confronted. Indeed, working with large 
animal groups usually implies that experiments are just 
observations, with the unavoidable difficulties to track 
trajectories and without much of a control parameter to 
vary Q. Similarly, experiments on living cells (during 
development or wound healing or within bacteria and 
amoeba colonies) often involve the presence of exter- 
nal (chemical) gradients, genetic factors, etc., which are 
hard to evaluate and known to have a possibly strong 
influence. There is hence a crucial need for model ex- 
periments using man-made objects with rather well un- 
derstood interactions. Swimmers Q (from chemically- 
powered nanorods to microscopic and macroscopic size 
mechanical devices) offer an interesting direction but the 
, intrinsically long-range nature of hydrodynamic interac- 
tions may appear as a unnecessary complication. In this 
context, vibrated, dry, inert, "granular" particles appear 
as an attractive case where much control can be exerted 
on the system, in the absence of long-range interactions 
or unwanted additional features, so that the onset of col- 
lective motion would then be a bona fide spontaneous 
symmetry breaking phenomenon. 

Various objects can be set in fairly regular motion on 
a flat surface when vibrated properly: Yamada, Hondou, 
and Sano were pioneers in demonstrating that an axisym- 
metric polar object vibrated between two plates can move 
quasi-ballistically the collective level, Kudrolli's 




FIG. 1: Collective motion of self-propelled disks (color on- 
line). Bottom left: sketch of our polar particles. Main panel: 
snapshot of an ordered regime observed in our flower-shape 
domain. The color scale reflects the local polar order from 
perfect alignment (red) to anti-alignment (blue). The intrin- 
sic polarity of the particles is indicated by the black arrows. 



group studied the behavior of polar rods [ll| and, more 
recently, of short snake-like chains [Hj , but was unable to 
observe genuine long-range orientational order, i.e. col- 
lective motion. A few other works have dealt with the col- 
lective properties of shaken elongated apolar particles (a 
realization of so-called "active nematics") 131, but there 



no net collective motion is expected anyway. 

Thus, to our knowledge, no well-controlled experiment 
has produced a fluctuating, collectively-moving ordered 
phase of the type frequently observed in simple numerical 
models. This may be just due to the scarcity of attempts, 
but recent results might provide a deeper reason: it was 
found that self-propelled particles with apolar (nematic) 
alignment interactions cannot give rise to polar order, i.e. 
to collective motion [l^ . (They may give rise, however, 
to nematic order.) The few experiments mentioned above 
all dealt with elongated objects ("self-propelled rods"), 
and fall into this class because of their shape. 

In this Letter, we report on experiments conducted on 



vibrated disks with a built-in polar asymmetry which en- 
ables them to move coherently (Fig. [T]). The isotropic 
shape of the particles and their rather specific inelas- 
tic collision properties prevent strong nematic alignment. 
Varying the amplitude of the vibration, we observe the 
onset of large-scale collective motion and the existence 
of giant number fluctuations with a scaling exponent in 
agreement with the predicted theoretical value. We dis- 
cuss the difficulties in characterizing collective motion in 
a finite domain and the possible key differences with the 
simple models usually considered at the theoretical level. 

Experiments with shaken granular particles are noto- 
riously susceptible _to systematic deviations from pure 
vertical vibration 
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We use a 110 mm thick trun- 
cated cone of expanded polystyren sandwiched between 
two nylon disks. The top disk (diameter 425 mm) is cov- 
ered by a glass plate on which lay the particles. The 
bottom one (diameter 100 mm) is mounted on the slider 
of a stiff square air-bearing (C40-03100-100254,IBSPE), 
which provides virtually friction-free vertical motion and 
submicron amplitude residual horizontal motion. The 
vertical alignment is controlled by set screws. The vibra- 
tion is produced with an electromagnetic servo-controlled 
shaker (V455/6-PA1000L,LDS), the accelerometcr for 
the control being fixed at the bottom of the top vi- 
brating disk, embedded in the expanded polystyren. A 
400 mm long brass rod couples the air-bearing slider 
and the shaker. It is flexible enough to compensate 
for the alignment mismatch, but stiff enough to en- 
sure mechanical coupling. The shaker rests on a thick 
wooden plate ballasted with 460 kg of lead bricks and 
isolated from the ground by rubber mats (MUSTshock 
100xl00xEP5,Musthanc). We have measured the me- 
chanical response of the whole setup and found no res- 
onances in the window 70 — 130 Hz. Here, we use a 
sinusoidal vibration of frequency / = 115 Hz and vary 
the relative acceleration to gravity F = 2TTaf^ /g. The 
vibration amplitude a at a peak acceleration of 1 g at 
this frequency is 25 ^m. Using a triaxial accelerometcr 
(356B18,PCB Electronics), we checked that the horizon- 
tal to vertical ratio is lower than 10^^ and that the spatial 
homogeneity of the vibration is better than 1%. 

Our polar particles are micro-machined copper- 
beryllium disks (diameter c? = 4 mm) with an off-center 
tip and a glued rubber skate located at diametrically op- 
posite positions (Fig. [Ij. These two "legs", which have 
different mechanical response under vibration, endow the 
particles with a polar axis which can be determined from 
above thanks to a black spot located on their top. Under 
proper vibration, they can be set in directed motion (see 
below). Of total height h = 2.0 mm, they are sandwiched 
between two thick glass plates separated by a gap of 
H = 2.4 mm. We also used, to perform "null case exper- 
iments" , plain rotationally-invariant disks (same metal, 
diameter, and height), hereafter called the "symmetric" 
particles. We laterally confined the particles in a flower- 
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FIG. 2: (color online) Individual dynamics for F = 2.7. (a) 
typical portions of polar particles trajectories inside the ROI. 
Black and red arrows indicate «' and n* at selected times. 
The domain area is about 15 x 15d. (b) same for symmetric 
particles, (c) pdf (lin-lin) of a, the angle between ii* and tI*. 
(d) variation of angular diffusion coefficient De with F. 

shaped arena of internal diameter D = 160 mm (Fig.[T]). 
The petals avoid the stagnation and accumulation of par- 
ticles along the boundaries as reported for instance in [llj 
by "reinjecting" them into the bulk. A CCD camera with 
a spatial resolution of 1728 x 1728 pixels and standard 
tracking software is used to capture the motion of the 
particles at a frame rate of 20 Hz. In the following, the 
unit of time is set to be the period of vibration and the 
unit length is the particle diameter. Within these units, 
the resolution on the position r of the particles is better 
than 0.1, that on the orientation n is of the order of 0.05 
rad and the lag separating two images is tq = 5.75. Mea- 
suring the long-time averaged spatial density map (for 
various numbers of particles), we find that this density 
field slightly increases near the boundaries, but is con- 
stant to a few percent in a region of interest (ROI) of 
diameter 2Qd. This provides an additional check of the 
spatial homogeneity of our setup. 

We first performed experiments with 50 particles, i.e. 
at a surface fraction small enough so that collisions 
are rare and the individual dynamics can be investi- 
gated. For large acceleration, the polar particles de- 
scribe random-walk like trajectories with short persis- 
tence length. Decreasing F, they show more and more 
directed motion, and the persistence length quickly ex- 
ceeds the system size. This is in contrast with the sym- 
metric particles which retain the same shortly correlated 
individual walk dynamics for all F values (Fig.[2^b). 

More precisely, individual velocities i7i(<) = {ri{t+To) — 
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FIG. 3: (color online) Trajectories of two particles "during" 
a collision: they first collide almost head-on, but repeated 
contacts (all along the red arrow) finally leave them almost 
aligned, despite their isotropic shape. 



fi{t))/To measured within the ROI have a well-defined 
most probable or mean value t^typ — 0.025 which changes 
by only 6% over the interval T G [2.7,3.7] (not shown). 
For smaller values of F the velocity decreases suddenly 
and the particles come to an almost complete stop around 
F = 2.4. The local displacements of our polar parti- 
cles are overwhelmingly taking place along ni(t), their 
instantaneous polarity (Fig. [2};). The distribution of the 
angle 9i(t,t + tq) by which they turn during an inter- 
val To (defined using the polarity ni{t)) is an exponential 
distribution of zero- mean and variance 2Dg/To. The an- 
gular diffusion constant Dg decreases fast and linearly for 
F e [2.7,3.7] (Fig. In contrast again, De is about 

one order of magnitude larger for our isotropic particles, 
and varies little with F (not shown) . A persistence length 
can then be defined as ^ = ^Tr'^vtyp/Dg (i.e. the length 
traveled over the time needed to turn by tt, assuming 
a constant speed ftyp)- Its typical value decreases from 
above 100 for F = 2.7 to around 20 for F ~ 3.7 whereas 
it stays around 1 for the symmetric particles. 

Wc now turn to the collective dynamics of our polar 
particles. As seen above, the relative acceleration F has 
a strong influence on their individual dynamics, control- 
ling the persistence length of their trajectories via the 
angular diffusion constant Dg. During collisions, they 
typically bounce against each other several times, yield- 
ing, on average, some degree of alignment (Fig. [J). All 
this is reminiscent of Vicsek-like models, for which one 
of the main control parameters is the strength of the 
angular noise competing with the alignment interaction 
0,|j]. Thus F is not only an easy control parameter, but 
also a natural one, which we use in the following. The 
surface fraction cj) of particles is another natural control 
parameter in collective motion and granular media stud- 
ies, but it is somewhat more tedious to vary, and, more 
importantly, one should avoid to deal with too few, re- 
spectively too many, particles in order to prevent loss 
of statistical quality, respectively jamming effects. Be- 
low, we present results obtained with N = 890 particles, 
which gives a surface fraction </> ~ 0.38 in the ROI where 
an average of 160 particles (slightly dependent on F) is 
found. Similar results were obtained at nearby densities. 
To characterize orientational order, we use the modulus 
of the average velocity-defined polarity 4'(t) ~ \{ui{t))\ 
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FIG. 4: (color online) Collective dynamics, (a) time series 
of order parameter * at F = 2.8. (b) pdf (lin-lin) of *(t) 
at various F values, (c) (*!') vs F for polar and symmetric 
particles, (d) An vs n at F = 2.8. 

where Ui (t) is the unit vector along Vi {t) and the average 
is over all particles inside the ROI at time t. [l^ 

At low F values, for which the directed motion of our 
polar particles is most persistent, we observe spectacu- 
lar large-scale collective motion, with jets and swirls as 
large as the system size (Fig. [T] and [13 )• Of course, 
because our boundary conditions are not periodic, the 
collective motion observed is not sustained at all times. 
Large moving clusters form, then breakdown, etc. As a 
result, the times series of the order parameter 'I' presents 
strong variations, but can take a rather well-defined order 
one value for long periods of time (Fig. 0^) . At high F 
values (large noise) no large-scale ordering is found. De- 
creasing F, the pdf of ^ becomes wider and wider, with 
a mean and a most probable value increasing sharply 
(Fig. IIIdc). Note that the most probable value corre- 
sponds, at small F, to the plateau value found in time 
series of \I'. Thus, we observe the clear emergence of 
long-range orientational order over the range of usable F 
values. In contrast, the same experiments realized with 
our symmetric particles do not give rise to any collective 
motion (Fig.|4j:), which ultimately indicates that our ob- 
servations with polar particles are not due to some resid- 
ual large-scale component of our shaking apparatus [15 1 . 

Unfortunately, we could not observe the saturation of 
the order parameter expected deep in the ordered phase, 
because the "self-propulsion" of our polar particles dete- 
riorates for F < 2.7. Nevertheless, large \E' values were ob- 
served, signalling that our lowest usable F values are al- 
ready in the ordered phase, albeit not quite surely out of 
the critical/transitional region. We thus investigate the 
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emergence of the so-called "giant number fluctuations" 
(GNF) which have been shown theoretically and numer- 
ically to be a landmark of orientationally-ordered phases 
for active particles [Hill- To this aim, we recorded, along 
time, the number n{t) of particles present in boxes of var- 
ious sizes located within the ROI. GNF are characterized 
by the fact that the variance An of this number scales 
faster than the mean n. This is indeed what we find: over 
a range of scales. An grows like n" with a ^ 1.45 ± 0.05. 
For larger scales, one feels the finite system size and An 
levels off. In fact, according to the prediction derived 
from the work of Tu and Toner [Tsj and confirmed in sim- 
ulations this number should be 1.6. Thus our finding 
is quite consistent with the predicted value, all the more 
so since a is expected to converge from below as the 
system size increases Q. Although this will require to 
be confirmed by experiments performed in larger dishes, 
this result constitutes the first experimental evidence for 
GNF in collections of polar active particles [l^ . 

To summarize, we have shown that shaken particles 
with a polarity not related to their shape can exhibit col- 
lective motion on scales or the order of the domain in 
which they evolve. In the most ordered regimes reach- 
able, we recorded giant number fluctuations with a scal- 
ing exponent consistent with that of polar active phases. 

That we observe dominant polar order is worth dis- 
cussing: it was recently shown that if their alignment 
interaction is nematic, polar particles cannot order po- 
larly and only nematic order arises, even if the particles 
are pointwi se |14| . This nematic order is made of po- 
lar packets [20[, which could dominate the global order 
in a small domain such as our ROI. It is not clear, at 
this point, whether our system falls in this class. As a 
matter of fact, we do observe a small fraction of parti- 
cles going "against" the main flow in our most ordered 
regimes (Fig. [l]and [l3|)- But this remains rare — most 
of the time, polar alignment is observed, see Fig. [3] — , in 
contrast to the above-mentioned studies, but in line with 
the numerical work of Grossman, Aranson, and Ben Ja- 
cob 
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Further investigations will require a detailed 
study of the statistics of collisions. 

In Vicsek-style models (and their continuous descrip- 
tions) no GNF proper exist near the transition, where 
high-order high-density bands emerge 0, [l^ . Here, we 
found GNF in our most-ordered regimes, with approxi- 
mately the expected exponent. But it is impossible, at 
this stage, to disentangle fluctuations due to the prox- 
imity of the transition, those due to the frustration in- 
duced by our boundaries (which would break bands) , and 
"genuine" GNF. Thus, performing experiments in larger 
domains and looking for parameter values which would 
allow to go deep into the ordered phase is of utmost im- 
portance. This could also allow to study the nature of 
the transition to collective motion. To this aim, an even 
better control of our vibration table is necessary, a task 
we are currently pursuing. 
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